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We compute gluon and quark propagators in 2-colour QCD at large baryon chemical potential jj.. 
The gluon propagator is found to be antiscreened in the superfluid, confined phase and screened 
in the large-/!, deconfined phase. We present the first attempt to determine corresponding elec- 
tric and magnetic gluon masses. The quark propagator undergoes dramatic modifications in the 
superfluid region as a result of the formation of a superfluid gap. These modifications include the 
appearance of zero crossings in the vector part of the (normal) quark propagator, a large suppres- 
sion of the scalar part, and the emergence of anomalous propagation. 
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1. Introduction 

Determining the phase diagram of QCD at large baryon density and small temperatures re- 
mains one of the outstanding problems of strong interaction physics. This problem is of both theo- 
retical and phenomenological interest: on the theoretical side, an exceptionally rich phase structure 
may be present, while the phenomenological interest is spurred by the possibility that some of these 
phases may be present in compact stars, and may have observable consequences. 

Direct lattice simulations of QCD at high density and low temperature are hindered by the sign 
problem, so alternative approaches are required. One such approach is to study QCD-like theories 
which may be simulated on the lattice, and apply the lessons learnt from these theories to the case 
of real QCD. Foremost among these theories is QCD with gauge group SU(2) (QC2D). 

Medium modifications of quark and gluon propagators is one topic where QC2D may directly 
inform real QCD calculations. The gluon propagator is used as input into the gap equation for the 
superfluid gap at high density, but the propagator that is used is usually based either on (resummed) 
perturbation theory or on simple generalisations of the vacuum propagator. Nontrivial medium 
modifications or nonperturbative effects may thus significantly alter the results. The quark propa- 
gator encodes information about effective quark masses and gap parameters, while first-principles 
results for gluon and quark propagators together can be used to check the assumptions going into 
dense QCD calculations in the Dyson-Schwinger equation framework [jlj 

2. Formulation 

We will be using Nf = 2 degenerate flavours of Wilson fermion, with a diquark source j 
included to lift low-lying eigenvalues and study diquark condensation without uncontrolled ap- 
proximations. The fermion action can be written 



where M(/l) is the usual Wilson fermion matrix with chemical potential ju. It satisfies the symme- 
tries 



with K = C75T2. The first of these is the Pauli-Giirsey symmetry. The inverse of M is the Gor'kov 
propagator, 




(2.1) 
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propagator has the symmetry properties 
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We will also write the inverse propagator as 




(2.6) 



which has the same symmetry properties as . 

The normal propagator S can in general be written in terms of four momentum-space form 
factors, 

S(p) =$S a (p 2 ,p 4 ) + S b (p 2 ,p 4 ) + Ja{pa ~ in)S c {p, P4 ) + ijA $S d (f,P4) ■ (2-7) 

In QC2D the Pauli-Gursey symmetry ensures that all form factors are purely real. The structure 
of the anomalous propagator depends on the pattern of diquark condensation. Assuming that the 
condensation occurs in the colour singlet channel with quarks of unequal flavour, the anomalous 
propagator can be written as T(p) = T'{p)CTT2 (and similarly for the anomalous part A(p) of the 
inverse propagator), where T = 75 for condensation in the scalar (0 + ) channel. Spin-1 condensation 
leads to more complicated structures, but is energetically disfavoured compared to spin-0 conden- 
sation and will not be considered here. The remaining spin structure can be written in terms of 
form factors T a ,Tb,T c ,T,i analogous to ([O]), ie 

T'ip) =pT a (p 2 , P4 ) + T b {p 2 ,p A ) + Ya{pa ~ ipl)T c (p,p 4 ) + ijA fT d (f,PA) ■ (2-8) 

Similarly, the inverse propagator can be written in terms of form factors A, fi, C and D for the normal 
part N, and <j> a ,<j>b, <t> c ,<t>d for the anomalous part A'(p). The form factors </>,- are the gap functions. 

The gluon propagator in presence of a chemical potential in Landau gauge may be decomposed 
into an magnetic and electric form factor, 

£> M v (q, qo) = P£ v D M (f ,ql)+ P^yD E (f , <&) ■ (2-9) 

The projectors P^ v (q),P^ v {q) are both 4-dimensionally transverse, and are spatially transverse and 
longitudinal respectively. 



3. Results 

We have generated gauge configurations on two lattices: a "coarse" lattice with j3 = 1.7, K = 
0.178, V = 8 3 x 16, and a "fine" lattice with /3 = 1.9, K = 0.168, V = 12 3 x 24. The lattice spacings 
are 0.23fm and 0. 18fm respectively, while m n jm p = 0.8 in both cases. A range of chemical poten- 
tials \i were used with diquark source a] = 0.04, while additional configurations were generated 
with a] = 0.02,0.06 for selected values of \i. In addition to this, we have also generated configura- 
tions at jl = for two "finer" lattices, with j8 = 2.0, fc = 0.162 ("heavy") and K = 0.163 ("light"), 
both with volumes V = 12 3 x 24. 

3.1 Gluons 

Results for the gluon propagator on the coarse lattice have been presented in 0; we will 
supplement those here with results from the fine lattice. On both lattices, an onset transition to a 
phase with nonzero baryon density and diquark condensate was found at /i « m K /2, while BCS- 
like scaling of energy density, baryon density and diquark condensate was found at higher jj.. On the 
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Figure 1: The gluon dressing function at zero chemical potential, for different lattice spacings and volumes. 



coarse lattice the crossover to BCS-like scaling was associated with a nonvanishing Polyakov loop 
L, indicating a coincident deconfinement transition [Q]. On the fine lattice, these two transitions are 
separate, with the deconfining transition occuring at considerably larger /I [Q]. 

First of all, we investigate the scaling behaviour of the gluon propagator in the vacuum (jj. = 0). 
Figure [j] shows the gluon dressing function q 2 D(q) for three of our four different lattices. For the 
coarse lattice parameters, we also have data for two different volumes. The data have all been 
cylinder cut to select the points with smallest lattice artefacts. Since the lattice spacing for the 
finer lattice has not yet been independently determined, the matching procedure described in [||] 
has been used to find the ratio of lattice spacings cif/a/f that gives the best match for the gluon 
propagator on the fine (/) and finer (ff) lattices. 

We see that finite volume effects are small for the momenta considered here, but scaling viola- 
tions (finite lattice spacing effects) are very large between the coarse and the two finer lattices. The 
good scaling observed between the two finer lattices may be somewhat misleading, since the match- 
ing procedure used in setting the scale for the finer lattice assumes we are in the scaling regime. 
Nonetheless, the good agreement over a wide range of momenta indicates that lattice artefacts here 
are not too large. 

Figure || shows the two lowest Matsubara modes of the unrenormalised gluon propagator as a 
function of spatial momentum \q\ for a range of chemical potentials, on both lattices. In all cases, 
the propagator at the lowest chemical potential /I shown is consistent with the vacuum propagator. 
On the coarse lattice both magnetic and electric propagator are strongly screened at large ju, while 
they are enhanced at low momentum in the intermediate-density region. The static (go = 0) mag- 
netic gluon propagator turns out to have a surprisingly strong dependence on the diquark source, 
which counteracts the infrared suppression at large /I as j — ► 0, but does not remove it completely. 
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Figure 2: The unrenormalised gluon propagator on the coarse lattice (top) and on the fine lattice (bottom), 
for various chemical potentials afi = 0.25 — 1.10. 



This is demonstrated in fig. [3L which shows the magnetic gluon propagator for the two lowest 
Matsubara frequencies, extrapolated to zero diquark source. We clearly see a strong infrared en- 
hancement at intermediate [l, but at a\i = 0.9(/J. = 0.78GeV) both the static and non-static modes 
are suppressed in the infrared. 

The same qualitative picture can be seen on the fine lattice, but in this case the infrared sup- 
pression sets in at much larger pL (around a\i = 0.8 or \i = 0.9GeV). This is consistent with the 
hypothesis that the screening effect is linked with the deconfinement transition, ie that it is a result 
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Figure 3: Magnetic gluon propagator on the coarse lattice, extrapolated to zero diquark source j. The 
left-hand plot shows the lowest Matsubara mode (q\ = 0), while the right-hand plot shows the first nonzero 
Matsubara mode. 



of the gluons being screened by coloured quark degrees of freedom. 

It is worth pointing out that the enhancement resp. screening noted here is in comparison to the 
vacuum gluon propagator, which is known to be infrared suppressed due to nonperturbative effects 
(as discussed at length in other contributions to this conference). It seems reasonable to assume 
that although the static magnetic gluon is unscreened to all orders in perturbation theory, nonper- 
turbative effects may be responsible for the additional screening observed here in the deconfined 
phase. 

We have attempted to fit the gluon propagator to a simple massive form, 

D EM (q,q 4 ;n) = Z ' ,m . (3.1) 

The resulting electric and magnetic gluon masses m e>m are shown as functions of pL in figure |j. It 
is worth noting that the quality of these fits is quite poor. This is expected, as it is known that at 
/i = it is not possible to describe the gluon propagator by a simple, momentum-independent mass, 
while at large /I one should reproduce the results of hard-dense-loop (HDL) resummed perturbation 
theory, which also has a more complicated functional form. A form which interpolates between 
HDL and available results for ji = ^, 0] is likely to yield better results. A further technical 
complication is that we have defined De only at nonzero spatial momenta, while the fits to Dm 
include the q = point. This is the reason for the discrepancy between m e and m m at /I = 0, where 
they should be equal. This also tends to yield lower values for m e throughout. 

With these provisos, we can see that both the electric and magnetic gluon masses remain 
roughly constant for small and intermediate jit, before increasing at large ju, corresponding roughly 
to the deconfined phase. We see, however, that there is a large difference between the mass values 
from the two lattices, indicating that scaling violations are still very large at these lattice spacings. 
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Figure 4: The electric and magnetic gluon mass as a function of chemical potential p., determined from a fit 
to a simple massive propagator on each lattice. For the coarse lattice, the filled symbols denotes fits to data 
with zero diquark source j, while the open symbols are from ja — 0.04. For the fine lattice all data are for 
ja = 0.04. It was not possible to get any fit for the electric gluon on the fine lattice for ajj, < 0.7. 



3.2 Quarks 

In the vacuum, there are only two independent tensor components of the quark propagator, 
which is conventionally written as 



Sip) 



(3.2) 



where M is the mass function and Z the renormalisation function. These are shown in fig. ^ 
for the different lattice spacings and quark masses available. Both Z(p) and M(p) have been 
multiplicatively tree-level corrected [J|]; however, since the critical quark mass is not yet known, 
the tree-level correction of M{p) is not yet properly carried out. 

We immediately see that there are large scaling violations in both form factors, and large 
violations of rotational symmetry in Z(p). In particular, we note that Z(p) increases in the infrared 
for the coarser lattices, whereas it is usually found to be infrared suppressed. We see that this 
suppression appears to be recovered as we move towards the continuum limit. A careful continuum 
extrapolation will be needed to obtain quantitative results. 

At nonzero chemical potential, we find that the form factors S a , Sb and S c (spatial- vector, scalar 
and temporal-vector) of the normal quark propagator and the form factors Tb and Tj (scalar and 
tensor) of the anomalous propagator are nonzero, while the remaining components are zero. Results 
for the coarse lattice were shown in [Q] ; here we will show results for the fine lattice only. 

Figure H shows the spatial- vector part S a and scalar part Sb of the normal quark propagator for 
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p [GeV] p [GeV] 

Figure 5: The quark propagator renormalisation function (left) and mass function (right) at zero chemical 
potential, for different lattice spacings. 

a range of chemical potentials a}i = 0.25 — 1.1. These both exhibit dramatic medium modifica- 
tions. The scalar propagator Sb is strongly suppressed in the superfiuid phase, suggesting a drastic 
reduction in the in-medium effective quark mass. This is linked to the appearence of the diquark 
condensate: the chiral condensate rotates into the diquark condensate in the superfiuid phase Jl0|]. 
We would therefore expect to find the missing strength in the anomalous propagator. The change 
in behaviour is sudden and takes place around \i Q = m n /2, while for larger /I there is little change. 

The spatial- vector propagator S a is also infrared suppressed at large }i, but this suppression 
happens gradually as a function of pL, and sets in only above \i . At the largest densities we see that 
S a (k,k4 = TiT) becomes negative for small spatial momentum |jfc|. 

The two lowest Matsubara modes of the temporal- vector propagator S c are shown in [7|. We see 
that the lowest Matsubara mode (£4 = nT) becomes negative at intermediate momenta, approaching 
zero from below at high momenta. This is a dramatic change compared to the vacuum propagator, 
which stays positive at all momenta, and indicates the formation of a superfiuid gap. The location 
of the zero crossing in the A4 — » limit corresponds to the Fermi momentum kp. In accordance 
with this, the zero crossing moves to larger \k\ as jU increases. On closer inspection, we find that the 
second Matsubara mode (&4 = 3nT) also becomes negative for large /I (at large spatial momentum). 
It would therefore in principle be possible to extrapolate this zero crossing to £4 = and thus find 
kp as a function of pL. 

Figure || shows the nonzero components of the anomalous Gor'kov propagator. The dominant 
part is, as expected, the scalar part 7/,, but a clear signal is also found for the tensor part T4, in 
accordance with what was found on the coarse lattice [Qj. We find that the lattice artefacts in 
the scalar part is substantially reduced compared to the coarse lattice, while the tensor part is 
still subject to very large violations of rotational symmetry. It may therefore be open to question 
whether this component will survive the continuum limit. 
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Figure 6: The lowest Matsubara frequency of the spatial-vector (left) and scalar (right) part of the normal 
quark propagator, on the fine lattice, for different chemical potentials ji. 
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Figure 7: The temporal-vector part of the normal quark propagator, on the fine lattice, for different chemical 
potentials fj,. 



The scalar anomalous propagator shows a clear change in behaviour as one goes from small to 
large chemical potential. Firstly, we note that it increases between a\l = 0.25 and 0.35. The former 
point is below the superfluid transition, but anomalous propagation is present due to the explicit 
diquark source. We expect that 7), (and Td) will vanish in the 7 — > limit for pL < As \i increases 
above Tj, develops a plateau at low momentum, which extends to larger \k\ with increasing }X. 
At large /I, 7], thus becomes approximately constant, suggesting that anomalous propagation may 
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Figure 8: The scalar (left) and tensor (right) part of the anomalous quark propagator, on the fine lattice. 



be described by a momentum-independent diquark gap A. 



4. Discussion and outlook 

We have found substantial modifications of both gluon and quark propagators in the dense 
medium. In the superfluid, confined phase, the electric and magnetic gluon propagators are both 
enhanced in the infrared compared to the vacuum. In the deconfined phase, they are both screened 
(infrared suppressed). This screening is evident even in the static magnetic gluon, which is un- 
screened to all orders in perturbation theory. If these results carry over to SU(3) they would inval- 
idate the use of an unscreened static-magnetic gluon propagator in the gap equation at large }i. A 
careful analysis at different volumes and lattice spacings is however necessary to draw quantitative 
conclusions. 

The dramatic modifications seen in the quark propagator are directly related to the appearance 
of a diquark gap. Our next step will be to compute the form factors, including the diquark gap 
and mass function, by inverting the quark propagator. Further quantitative studies will include 
determining the Fermi momentum pp by extrapolating the zero crossing in the temporal-vector 
propagator S c to = 0, and determining the size of Cooper pairs from the anomalous propagator, 
to study the BEC-BCS crossover in more detail. 
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